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The intensity of Galactic cosmic rays is nearly isotropic due to the influence 
of magnetic fields in the Milky Way. Here we present two-dimensional high- 
precision anisotropy measurement for energies from a few to several hundred 
TeV using the huge data sample of the Tibet Air Shower Arrays. Besides re- 
vealing finer details of the known anisotropies, a new component of sidereal 
time Galactic cosmic ray anisotropy is uncovered around the Cygnus region 
direction. For cosmic-ray energies up to a few hundred TeV, all components of 
anisotropies fade away, showing a corotation of Galactic cosmic rays with the 
local Galactic magnetic environment. These results bear broad implications 
to cosmic rays, supernovae, magnetic field, heliospheric and Galactic dynamic 
environment in a comprehensive manner. 
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Galactic cosmic rays 

The anisotropy of Galactic cosmic rays (GCRs) may result from an uneven distribution of cos- 
mic ray (CR) sources and the process of CR propagation in the Milky Way. CRs of energy 
below 10 15 eV are accelerated by diffusive magnetohydrodynamic (MHD) shocks (0|2]|5]H1) of 
supernova remnants (SNRs) and stellar winds. The discreteness of SNRs could lead to a CR 
anisotropy $5^. However, GCRs must almost completely lose their original directional infor- 
mation; their orbits are deflected by the Galactic magnetic field (GMF) and are randomized by 
irregular GMF components, having traveled on average for many millions of years, some also 
having interacted with interstellar gas atoms and dusts. The transport of CRs in a magnetized 
plasma is governed by four major processes: convection, drift, anisotropic diffusion and adi- 
abatic energy change (deceleration or acceleration) (00. High-precision measurement of the 
CR anisotropy provides a means to explore magnetic field structures and gains insight for the 
CR transport parameters ©. The long-term high-altitude observation at the Tibet Air Shower 
Arrays (referred to as the Tibet AS7 experiment) has accumulated tens of billions of CR events 
in the multi-TeV energy range, ready for an unprecedented high-precision measurement of the 
CR anisotropy as well as the temporal and energy dependence of the CR anisotropy. 

An expected anisotropy is caused by the relative motion between the observer and the CR 
plasma, known as the Compton-Getting (CG) effect (9) with CRs arriving more intense from 
the motion direction and less intense from the opposite direction. Such a CR anisotropy, caused 
by the Earth orbital motion around the Sun, has indeed been detected HOWllh . Data assembled 
to 1930s © were consistent with a scenario that the CR plasma stays at rest in an inertial frame 
of reference attached to the Galactic center. If this were true, the Galactic rotation in the solar 
neighborhood might then be measurable. Nevertheless, such CR anisotropy due to the solar 
system rotation around the Galactic center at a speed of ~ 220 km s -1 remained inconclusive 
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over seven decades. Now our high-precision two-dimensional (2D) measurement gives a strong 
evidence to exclude the CR anisotropy of this origin and thus show a corotation of GCRs with 
the local GMF environment. 

Historically, the GCR anisotropy (f721[7?t has been measured as the sidereal time varia- 
tion at the spinning Earth using both underground fi detectors and ground-based air shower 
arrays d7?ll75ll7o"lli71l7ffll7Pl) . Located at different geographical latitudes and operating in dif- 
ferent years with various threshold energies, each individual experiment could only measure 
the CR modulation profile along the R.A. direction which was usually fitted by first few har- 
monics. Instead of using sine or cosine harmonics, one may adopt two Gaussian functions with 
Declination dependent parameters (mean, width and amplitude) to fit the so-called "tail-in" and 
"loss-cone" features, and a tentative 2D anisotropic picture was obtained (1201127 1) by simultane- 
ously fitting different experimental data. Here, the CR deficiency was thought to be associated 
with a magnetic cone-like structure and thus the name 'loss-cone', while the CR enhancement 
is roughly in the direction of the heliospheric magnetotail and is thus referred to as 'tail-in' 
enhancement (\12H13l . However, the accurate spatial and energy dependence of CR anisotropy 
could not be given (|22b . and some subtle features remain hard to reveal, since the CR anisotropy 
is more complex and cannot be properly described by two Gaussians. The Tibet AS-7 experi- 
ment alone can achieve 2D measurement in various energy ranges and to provide details of the 
2D CR anisotropy. 

Tibet Air Shower Array Experiment 

The Tibet Air Shower Array experiment has been conducted at Yangbajing (90.522 E, 30.102 
N; 4300 m above the sea level) in Tibet, China since 1990. The Tibet I array (l2~7t. consisting of 
49 scintillation counters and forming a 7 x7 matrix of 15 m span, was expanded to become the 
Tibet II array with an area of 36,900 m 2 by increasing the number of counters in 1994. In 1996, 
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part of Tibet II with an area of 5175 m 2 was upgraded to a high-density (HD) array with a 7.5 
m span (|2?b . To increase the event rate, the HD array was enlarged in 1999 to cover the central 
part of Tibet II as Tibet III array (1251 1261 \Ffa . The area of Tibet III array has reached 22,050 
m 2 . The trigger rates are ~ 105Hz and ~ 680Hz for the Tibet HD and III arrays, respectively. 
The data were acquired by running the HD array for 555.9 live days from 1997.2 to 1999.9 and 
the Tibet III array for 1318.9 live days from 1999.11 to 2005.10. GCR events are selected, if 
any four- fold coincidence occurs in the counters with each recording more than 0.8 particles 
in charge, the air shower core position is located in the array and the zenith angle of arrival 
direction is < 40°. With all those criteria, both Tibet HD and III arrays have the modal energy 
of 3 TeV and a moderate energy resolution; the ~ 0.9° angular resolution estimated from Monte 
Carlo simulations was verified by the Moon shadow measurement (|25]|26][27|). In total, 

~ 37 billion CR events are used in our data analysis. 

Data analysis and results 

With such a large data sample, we conduct a 2D measurement to reveal detailed structural 
information of the large-scale GCR anisotropy beyond the simple R.A. profiles. For each short 
time step (e.g. 2 minutes), the relative CR intensity at points in each zenith angle belt can be 
compared and this comparison can be extended step by step to all points in the surveyed sky [see 
Refs. (30) for details of data analysis]. Lacking the absolute detector efficiency calibration in the 
Dec direction, absolute CR intensities along different Dec directions cannot be compared. Thus, 
the average intensity in each narrow Dec belt is normalized to unity. Our analysis procedure 
would give a correct 2D anisotropy if there is no variation in the average CR intensity for 
different Dec. We systematically examined the CR anisotropy in four different time frames, 
namely solar time for solar modulation, sidereal time for Galactic modulation, anti- sidereal 
time and extended-sidereal time; and systematic variations are found to be unimportant. 
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To study the temporal variation of CR modulation, we divide the data sample into two sub- 
sets. The first subset is from 1997.2 to 2001.10 covering the 23rd solar maximum (a period 
of a few years when solar magnetic activity are the strongest) while the second subset is from 
2001.12 to 2005.1 1 approaching the solar minimum (a period of a few years when solar mag- 
netic activity become minimal). Comparing the sidereal time plots for these two intervals (Fig. 
1) shows that the CR anisotropy is fairly stable and insensitive to solar activity. The 'tail-in' and 
'loss-cone' anisotropy components (\12H13l . extracted earlier from a combination of the under- 
ground fi telescope data analyses (L20H2il) . are seen in our 2D plots in much finer details and with 
a high significance (Fig. 1). Our new high-precision measurement thus provides constraints on 
physical interpretations of these features. 

Spreading across ~ 280° to ~ 360° in R.A., a new excess component with a ~ 0.1% in- 
crease of the CR intensity peaked at Dec ~ 38° N and R.A. ~ 309° in the Cygnus region is de- 
tected at a high significance level of 13.3cr with a 5° pixel radius (Fig. Id). The Cygnus region, 
where complex features are revealed in broad wavelength bands of radio, infrared, X-ray and 
gamma-rays, is rich of candidate GCR sources. Recently, the first unidentified TeV gamma-ray 
source was discovered here by HEGRA (31). This region, as observed by EGRET (lJ2l . ap- 
pears to be the brightest source of diffuse GeV gamma-rays in the northern sky and contributes 
significantly to the diffuse TeV gamma-ray emission in the Galactic plane as observed by Mi- 
lagro (33) which rejects 90% of CR background while still retains ~45% of gamma-rays. Such 
gamma-rays originate from the interaction of CRs with gas and dusts. Using more stringent 
event selection criteria (07J) . a deeper view of Cygnus region with a 0.9° pixel radius shows that 
the large-scale excess consists of a few spatially separated excesses of smaller scales superposed 
onto a large-scale anisotropy (Fig. le); these small-scale (~ 2°) excess favors the interpreta- 
tion that the extended gamma-ray emission from the Cygnus region contributes significantly to 
the overall excess in the region (0?t. As our experiment cannot yet distinguish gamma-rays 
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from the charged CR background, we cannot tell how much of this excess is to be attributed to 
gamma-rays and how much, if any, is associated with charged CRs d53T> . Such a determination 
requires upgrading the Tibet arrays for CR and gamma-ray discrimination. 

The solar time CR modulation was also stable (Fig. 2). We found that including events with 
fewer than 8 detector coincidences (lower energy events) resulted in much larger modulation 
amplitudes than those obtained when these events were excluded (higher energy events). To 
avoid this, high multiplicity events with coincident detector numbers > 8 were adopted (Fig. 
2). The observed dipole anisotropy agrees very well with the expected CG effect due to the 
Earth orbital motion around the Sun. Thus, heliospheric magnetic field and solar activity does 
not influence the multi-TeV CR anisotropy. 

Because of the stable nature of the sidereal time modulation, data from different years were 
combined to examine the energy dependence of CR anisotropy. Fig. 3 shows the variation of 
anisotropy for five groups of events according to their different primary energies. For primary 
energies below 12 TeV, the anisotropics show little dependence on energy while above this 
energy, anisotropics fade away, consistent with a CR isotropy of KASCADE (17) in the energy 
range of 0.7 — 6PeV. Contrary to the earlier suggestion (13), the 'tail-in' component remains 
still visible above 50 TeV in smaller regions. Since the multi-TeV GCRs, whose gyro-radii are 
hundreds or thousands of AUs, are not affected by the heliospheric magnetic field, it is clear 
that the GMF must be responsible for both 'tail-in' and 'loss-cone' modulations. 

As a result of a diminishing GCR anisotropy at high energies, we can test the CG anisotropy 
caused by the orbital motion of the solar system in our Galaxy, which would peak at (a = 
315°, 5 = 49°) and minimize at (a = 135°, 5 = -49°) with an amplitude of 0.35%. This would 
be a salient signal in a real 2D measurement. However, as explained earlier, the modulation 
along the Dec direction is partly lost. After applying the normalization procedure along each 
Dec belt, the expected CG anisotropy is distorted and apparently peaks at (a = 315°, 5 = 0°) 
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and forms a trough around (a = 135°, 5 = 0°) with a smaller amplitude of ~0.23% (Fig. 4). To 
avoid any contamination from the non-vanishing 'tail-in' and 'loss-cone' anisotropies (\12H13h 
when the primary energy is of ~ 300 TeV, the upper half of the surveyed CR intensity map (with 
Dec > 25°) is used to compare with the predicted Galactic CG effect of an amplitude ~0.16%. 
The fitted anisotropy amplitude is 0.03%±0.03%, consistent with an isotropic CR intensity. 
Therefore our observation exclude the existence of the Galactic CG effect with a high-degree 
confidence, assuming the absence of other cancelling effects. The null result of the Galactic CG 
effect implies that GCRs corotate with the local GMF environment. 

Discussion 

The observation of GCR anisotropy and diffuse gamma-ray emission plays an important role in 
probing sources and propagation of CRs. The detection of the new large-scale GCR anisotropy 
component and the indication of extended gamma-ray emission from the same mysterious 
Cygnus region allow us to connect the GCR acceleration site and propagation. A precision 
spectral and morphological determination of the extended gamma-ray emission would be our 
next pursuit. The existence of large-scale GCR anisotropies up to a few tens of TeV indicates 
that they are not related to the helio spheric magnetic field. It is conceivable that GMF has 
large-scale structures in the heliospheric neighborhood. 

As in many spiral galaxies, our Milky Way has large-scale differential rotations in stellar 
and magnetized gas disks with a GMF of a few micro Gauss. The GMF, GCRs and thermal gas 
have similar energy densities of ~ leV cm -3 and interact with each other dynamically. The 
corotation of the GCR plasma with the local GMF environment around the Galactic center is 
enforced by the Lorentz force as GCRs randomly scatter and drift in irregular GMF compo- 
nents (OoT) . As the Galactic disk rotates differentially, the important inference is that the bulk 
GCR plasma within and above the Galactic disk must also rotate differentially. The GCR coro- 
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tation evidence provides an important empirical basis for the study of galactic MHD processes, 
such as modeling synchrotron emission diagnostics for large-scale spiral structures of MHD 
density waves QTTflMM- 
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Figure 1: Celestial CR intensity map* for Tibet AS7 data taken from 1997-2001 (a) and 2001-2005 (b). The 
vertical color bin width is 2.5 x 10~ 4 for the relative intensity in both (a) and (b). The circled regions labeled by 
I, II and III are the 'tail-in', 'loss-cone' (I72II73> and newly found anisotropy component around the Cygnus region 
(~ 38° N Dec and ~ 309° R.A.), respectively. Panel (c) is the one-dimensional (ID) projection of the 2D maps 
in R.A. for comparison. Panels (d) and (e) show the significance maps of Cygnus region (pixels in radius of 0.9° 
and sampled over a square grid of side width 0.25° for (e)) for data of 1997—2005. The vertical color bin widths 
are 0.69er and 0.42er for the significance in (d) and (e), respectively. Two thin curves in both (d) and (e) stand for 
the Galactic parallel b = ±5°. Small-scale anisotropics (e) superposed onto the large-scale anisotropy hint at the 
extended gamma-ray emission. ^Unless otherwise stated, images in Figs 1—4 are presented using pixels in radius 
of 5° and sampled over a square grid of side width 2°; the modal energy is 3 TeV. 
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Figure 2: Local solar time CR intensity map for the Tibet AS7 data taken from 1997-2001 (a) and (b) and 
2001—2005 (c) and (d). Both samples have the modal energy of 10 TeV. The vertical color bin width is 1.6 x 10~ 4 
for the relative intensity in both (a) and (c). In both (b) and (d), the fitting function is in the form of Amp * 
cos[27r(T— 4>)/2A] where the local solar time T and <f> are in unit of hour and Amp is the amplitude. The \ 2 
fit involves the number of degree of freedom (ndf) given by the number of bins minus 2 due to the two fitting 
parameters Amp and <f>. The ID plots are the projection of the 2D maps in local solar time. In the ID plots the 
dashed lines are from the expected CG effect, while the solid lines are the best harmonic fits, which agree very 
well with the prediction. The solar time modulation appears quite stable and insensitive to solar activity. 
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Figure 3: Celestial CR intensity map for different representative CR energies: (a) 4 TeV; (b) 6.2 TeV; (c) 12 TeV; 
(d) 50 TeV; (e) 300TeV. Data were taken during 1997-2005. The vertical color bin width is 2.5 x 10 -4 in (a)-(d), 
while it is 7.25 x 10~ 4 in (e) for different statistics, all for the relative CR intensity. 
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Figure 4: Celestial or 2D local sidereal time CR intensity map and its ID projection in the R.A. direction for 
300 TeV CRs of all data, (a) The colored map is the same as Fig. 3e, while the contours are the "apparent" 2D 
anisotropy expected from the Galactic CG effect. The vertical color bin width is 7.25 x 10 -4 for the relative 
intensity in (a). The ID projection is in map (b) for Dec between 25° — 70°, where the dashed line is the expected 
Galactic CG response while the solid line is the best fit to this observation using a first-order harmonic function. 
The fitting function is in the form of Amp * cos(R.A. — cf>) where <j> is in degree and Amp is the amplitude. The \ 2 
fit involves the ndf given by the number of bins minus 2 for the two fitting parameters Amp and <p. The data shows 
no Galactic CG effect with a high-level of confidence. 
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